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Abstract: We study one and two- loop triangle integrals with massless propagators and 
all external legs off shell. We show that there is a kinematic region where the results can 
be expressed in terms of a basis of single-valued polylogarithms in one complex variable. 
The relevant space of single-valued functions can be determined a priori and the results 
take strikingly a simple and compact form when written in terms of this basis. We study 
the properties of the basis functions and illustrate how one can easily analytically continue 
our results to all kinematic regions where the external masses have the same sign. 

Keywords: Feynman integrals, multi-loop computations, vertex diagrams. 



1. Introduction 



The evaluation of Feynman integrals is a necessary ingredient to higher-order corrections 
to physical observables in quantum field theories. For this reason analytic computations of 
two-loop Feynman integrals have seen a lot of interest over the last decade. The general 
strategy currently consists in reducing all the Feynman diagrams that contribute to a 
given process to small set of so-called master integrals by using integration-by-parts and 
Lorentz invariance identities [1, 2, 3, 4]. Various techniques have been developed over 
the last decade for the computation of the master integrals, the most prominent ones 
probably being the differential equation [5, 6, 7] and Mellin-Barnes approaches [8]. All 
of these procedures finally lead to analytical results for the master integral in the form 
of a Laurent series in the dimensional regulator e = (4 — D)/2 whose coefficients are 
transcendental functions of the scales of the process. While there are examples of multi- 
loop integrals which evaluate to elliptic integrals [9], it is known that for large classes 
of Feynman integrals these transcendental functions consist only in polylogarithms and 
generalizations thereof [10, 11, 12, 13, 14, 15, 16]. 

Focusing on two-loop 2 — >■ 2 processes relevant to hadron collider physics, all the 
two-loop master integrals with massless propagators for the production of two massless 
particles [17, 18, 19, 20, 21] or one massive and one massless particle [11, 22] have been 
computed analytically. Recently, master integrals for the production of a pair of heavy 
quarks at hadron colliders have also become available [23, 24, 25]. 

The focus of this paper are two-loop three-mass triangle integrals without internal 
masses. These integrals are for example relevant for the computation of the two-loop 
corrections to the production of two heavy particles with different masses. A prominent 
example of such a process is the production of a pair of weak gauge bosons at a hadron 
collider, an important background to many BSM and Higgs searches at the Large Hadron 
Collider (LHC). In ref. [26] it was shown that all two- loop three- mass triangle integrals with 
massless propagators can be reduced to a small set of master integrals consisting of one and 
two-loop two and three-point functions. In this paper we present analytical expressions 
for all the three-point master integrals presented in ref. [26]. While this is not the first 
time that these integrals have been evaluated analytically, we believe that our results go 
beyond existing representations available in the literature. Indeed, while some of the master 
integrals were already evaluated almost two decades ago in refs. [27, 28, 29] in terms of 
classical polylogarithms, not all of them were computed to the order in the e expansion 
required for two-loop computations in general, i.e., up to terms in the e expansion of 
transcendental weight four. In ref. [26] the two-loop master integrals were computed up to 
weight four using the differential equations approach and the results were expressed in terms 
of complicated iterated integrals whose integration kernels involve inverse square roots of 
a Kallen function. The resulting set of functions is then not related to polylogarithms in a 
straightforward way. 

The aim of this paper is to present fully analytic results for all two-loop three-mass 
triangle integrals without internal masses in terms of multiple polylogarithms up to weight 
four. We show that, using generic considerations on the analytic structure of these integrals, 
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one can identify a priori a basis for a space of poly logarithmic functions through which all 
of these integrals can be expressed. This space of functions is composed of single-valued 
functions of a single complex variable z (and its complex conjugate z), and encompasses 
in particular the famous Bloch-Wigner dilogarithm and the single-valued versions of the 
harmonic polylogarithms introduced in ref. [30]. The latter functions were recently used 
to derive analytic results for the six-point amplitude in N = 4 Super Yang-Mills in the 
multi-Regge limit [31] and for certain generalized ladder integrals [32] for a high number of 
loops. In terms of this basis of functions, all our results are characterized by very compact 
analytic expressions that make all the symmetries and analytic continuations manifest. 
In particular, we observe that one of the two-loop master integrals can be expressed, at 
least up to transcendental weight four, as a combination of the one-loop triangle and the 
two-loop ladder integral. This extends an observation made in ref. [28] to one order higher 
in the e expansion. 

This paper is organized as follows: In section 2 we define our notations and conventions 
and discuss some general properties of three-mass triangle integrals. In particular we argue 
that, in a specific kinematic region, they can naturally be expressed through a certain class 
of single- valued functions of a single complex variable z and its complex conjugate z, and 
in section 3 we give a short review of these functions up to weight four. In section 4 we 
present our results in the kinematic region where the functions are single- valued, and in 
section 5 we perform the analytic continuation of our results to other kinematic regions. 
We include appendices that contain details about the basis of single-valued functions. 



2. Triangle integrals with three external masses 

We start by discussing the kinematics of three-point functions where all three external legs 
are off shell and all internal propagator are massless. If T^(pi,P2 5 2>3; e) denotes a generic 
£-loop integral of this type in D = 4 — 2e dimensions with external momenta pi, i = 1,2, 3, 
then Lorentz invariance and momentum conservation imply that the result can only depend 
on the virtualities pf ^ 0. In dimensional regularization we can therefore write, without 
loss of generality, 

T^(p 1 ,p 2 , P3 ;e) = 4(-plr-^r^(u,v;e) , (2.1) 
for some integer n and where we defined 

2 2 

u = % and v = %. (2.2) 
Pi Pi 

In eq. (2.1) we pulled out the usual loop factor 

cr-e r(l-2e) ' {6) 

where 'Je = — r'(l) is the Euler-Mascheroni constant. In the following, and unless stated 
otherwise, we will always work in the Euclidean region where pf < 0, i = 1, 2, 3, and thus 
u, v > 0. All the results we present in this paper are real in the Euclidean region. We 
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Figure 1: The master integrals for the two- loop three- mass triangle integrals. 

note that our results will be equally valid in the physical region pf > 0. This region is 
phenomenologically relevant for two reasons. Firstly, it describes the decay of a heavy 
particle of mass p\ into two lighter particles 1 . Secondly, three-mass triangle integrals in 
this region appear in the production amplitude for a pair of weak gauge bosons at higher 
orders in perturbation theory. The region pf > is related to the Euclidean region via the 
analytic continuation 

_(p2 +ie )_> e -*r|p2| . (2.4) 

It is then easy to see that the phase factors cancel out in the ratios (2.2) so that the analytic 
continuation of eq. (2.1) from the Euclidean to the physical region p 2 > is trivial, 

T^(puP2,p 3 ;e) -)• {-l) n e*" u 4 \pir £t T®(u,v;e). (2.5) 

The main focus of this paper are three-mass triangles that appear in two-loop compu- 
tations in dimensional regularization. In ref. [26] it was shown that all two-loop three-mass 
triangles can be reduced to a limited set of master integrals (see fig. 1). The reduction to 
master integrals involves also two-loop two-point functions which can be evaluated to all 
orders in e in terms of T functions. We will therefore not consider the two-point master 
integrals any further and will concentrate exclusively on the three-point functions. 

As we will see in the next sections, the kinematics of a genuine three-point function is 
most conveniently parametrized in terms of two variables z and z, 

zz = u and (1 — z) (1 — z) = v , (2-6) 

or equivalently 

z=-(l + u — v+ \J X(l,u,v)^j and z=-(l + u — v— \j \(l,u,v)^j , (2.7) 

where A denotes the Kailen function 

X(a, b, c) = a 2 + b 2 + c 2 - 2ab - 2ac - 26c . (2.8) 

The appearance of the Kailen function in eq. (2.7) divides the (u, v) plane into four different 
kinematical regions, shown in fig. 2. In the regions II, III and IV the Kailen function 

1 Without loss of generality we may assume that p% is the largest invariant. 
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Figure 2: Different regions in (u,v) space. The parabola represents the curve X(l,u,v) = 0. 

X(l,u,v) is positive, and so z and z are both real. In region I, on the contrary, we have 
A(l, u, v) < 0, and thus z and z are complex conjugate to each other. There is a fifth region 
where the Kallen function vanishes, i.e., the boundary between region I and the regions 
II, II and IV. If we consider decay kinematics, region I is kinematically not allowed and 
has therefore often been discarded in the literature. In the following we take a different 
viewpoint and we argue that region I is the fundamental domain in which all three-mass 
triangle integrals are defined. The other regions are related to this fundamental domain 
by analytic continuation. In order to motivate this statement we need to introduce a 
mathematical tool that allows us to analyze the structure of polylogarithmic functions, the 
so-called symbol, which we briefly review in the rest of this section. 

2.1 Symbols of three-mass triangles and single-valued multiple poly logarithms 

Before defining the symbol map, we first need to define multiple polylogarithms, which are 
the natural set of functions describing large classes of Feynman integrals. Multiple poly- 
logarithms are a generalization of the ordinary logarithm and the classical polylogarithms, 



with Lii(z) = — ln(l — z). It is known that multi-loop multi-scale integrals can give rise 
to new classes of functions, among which one finds the so-called multiple polylogarithms, 
a multi-variable extension of eq. (2.9) defined recursively via the iterated integral [12, 13] 




(2.9) 




(2.10) 
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with G{z) = 1 and where a,i,z € C. In the special case where all the a^s are zero, we 
define, using the obvious vector notation a n = (a, . . . , a), 



G(O n ;z) = - } ln n z. (2.11) 
n! 

The number n of elements a,, counted with multiplicities, is called the weight of the multiple 
polylogarithm. In special cases multiple polylogarithms can be expressed through classical 
polylogarithms and ordinary logarithms only, e.g., 

G(a n ; z) = -. In™ (l - -) and G(O n _i, a; z) = -Li n (-) . (2.12) 
n! V a/ \a/ 

In general, however, no such simple formulae are known. An important role is played in 
physics by the so-called harmonic polylogarithms (HPLs) [10], which correspond, up to a 
sign, to the special case a* € {—1, 0, 1}, 

H(ai,...,a n ;z) = (-l) p G{a u . . . , a n ; z) , (2.13) 

where p denotes the number of elements in (a±, . . . , a n ) equal to +1. 

Iterated integrals form a shuffle algebra [33], which allows one to express the product 
of two multiple polylogarithms of weight n\ and ni as a linear combination with integer 
coefficients of multiple polylogarithms of weight n\ +ri2, 

G(ai, . . . , a ni ; z) G(a ni+ i, . . . , a ni+n2 ; z) = G(a a ^, . . . ,a a ^ ll+n2 y, z), (2 14) 

o-eS(ni,ri 2 ) 

where £(711,712) denotes the set of all shuffles of n\ + 712 elements, i.e., the subset of the 
symmetric group S ril+ri2 defined by 

£(711,712) = {cr e S ni+n2 \ a^ 1 (I) < ... < cr _1 (ni) and a^ 1 (n 1 + l) < ... < cr _1 (7ii+n 2 )} . 

(2.15) 

Multiple polylogarithms satisfy various complicated functional equations among them- 
selves. A way to deal with the functional equations is given by the symbol, a linear map 
which associates to every multiple polylogarithm an element in the tensor algebra over the 
group of rational functions. Various (equivalent) definitions have been given in the litera- 
ture for the symbol of a multiple polylogarithm [34, 35, 36, 37, 38, 39]. One possible way 
to define the symbol of a multiple polylogarithm is to consider its total differential [13], 



n-1 



n— 1 / 

dG(a n -i, . . .,ai;a n ) = ^G(a„_i, . . . ,a i: . . . ,ai;a n )dln ( 

i=i ^ 

and to define the symbol recursively by [37] 



O-i — Oi-\ 



(2.16) 



n— 1 / 

S(G(a n - 1 ,...,a 1 -a n )) = ^S(G(a n -i,...,a i ,...,a 1 ;a n ))® I l _ 1+1 ) . (2.17) 



i=l 
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As an example, the symbols of the classical polylogarithms and the ordinary logarithms 
are given by 

<S(Li n (z)) = -(l-z)8 z8. . .®z and S f ln n z] = z8. . .8z . (2.18) 

(n— 1) times n times 

In addition the symbol satisfies the following identities, 

...®(a-b)®... = ...<8>a<8>... + ...<8>&<8>... , 

2.19 

...®(±1)®... =0. 

The symbol also encodes information about the cuts of a function. As an example, let 
us take a function f(z) whose symbol can be written schematically in the form 

S(f(z)) = ai (z)®...®a n (z), (2.20) 

where a>i{z) are rational functions of some variable z. Then f{z) has branch cuts in the 
complex z plane starting at points Zi with a\(zi) = or a±(zi) = oo, and the symbol of the 
discontinuity across the branch cut is obtained by dropping the first entry in the symbol 
of f(z) , 

S [Disc ai(2 . )=o /0z)] = a 2 (z) ® . . . ® a n (z) . (2.21) 

If f(z) is a loop integral, then its cuts are determined by Cutkosky's rules. This implies in 
particular that the first entry of the symbol of a loop integral (with massless propagators) 
must be a Mandelstam invariant [40]. In the special case of three-mass triangles with 
massless propagators we consider, this implies that at each order in the e expansion, 

oo 

T (e) (u, v;e) = Y, 7~ fc W (u, v) e k , for some m € Z , (2.22) 

k=m 

the symbol of the coefficients of the Laurent series must take the form 

S rP{u,v)^ = ^2lR'?(u,v)u®Ui{u,v) + R"(u,v)v ®Vi{u,v)} (2.23) 

i 

= [Ri( z , *) ® Ui(z, z) + %{z, z) [(1 -z)(l- z)] ® Vi(z, z)] , 

i 

where R™' v are algebraic functions of the kinematic variables u and v, and U{ and Vi are 
tensors of lower weight. 

The form of the symbol (2.23) exhibits an important property of the functions Trf 1 ■ 
Let us assume that we are working in region I where A(l,u, v) < 0, and thus z and z 
are complex conjugate to each other. If we want to compute the discontinuities of the 
function T^p (seen as a function of (z,z)) in the complex z plane, then it follows from 
eq. (2.21) that Tfp has potential branch cuts starting at z = 0, z = 1 and z = oo. We 
will concentrate for now on the discontinuity around z = 0. The argument for the other 
two cases is similar. The symbol of Disc2=o7^, only has contributions from the first term 
in eq. (2.23). Terms of the form z <g> U{ and z ®Ui obviously contribute with opposite 
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signs, and so they cancel. The same argument holds for the discontinuities around z = 1 

(£) 

and z = oo. We thus conclude that the functions 7^ , seen as functions of the single 
complex variable z and its complex conjugate z, have no branch cuts in the complex z 
plane, and hence they are a single-valued functions of the complex variable z. We have 
thus shown that three-mass triangles without internal masses are expressible in region 
I through single-valued polylogarithms of the complex variable z. This is in agreement 
with the corresponding result for the generalized ladder integrals [32], which reduce to 
three-mass triangle integrals upon using conformal invariance to send a point to infinity. 
We note however that while, strictly speaking, single-valuedness implies the first entry 
condition (2.23), the inverse is not necessarily true. We discuss this problem in more detail 
in appendix A. It is however possible to carry out a more careful analysis based on the 
Hopf algebra of multiple polylogarithms without altering the conclusion. The requirement 
that the functions be single-valued then strongly constrains the set of functions that can 
appear as coefficients in the Laurent series (2.22). This class of functions will be reviewed 
in the next section. 

3. Single-valued multiple polylogarithms in one variable 

In this section we give a short review of single-valued polylogarithms in one complex vari- 
able z. We start by defining the single-valued analogues of the classical polylogarithms. 
The classical polylogarithms have branch cuts starting at z = 1, and the discontinuity 
across the branch cut is given by 



The knowledge of the discontinuities (3.1) can be used to construct linear combinations of 
classical polylogarithms in z and its complex conjugate z such that all the discontinuities 
cancel. As a result, we obtain a sequence of real- analytic functions on the punctured 
complex plane C/{0, 1}. Although the space of single- valued functions is unique, there 
is a freedom of how to choose a basis in the space of single- valued functions, and several 
definitions of single- valued analogues of the classical polylogarithms have been proposed in 
the literature. Here we use the definition of Zagier, and we define 



Kk=o ) 

where $\ n denotes the real part if n is odd and the imaginary otherwise, and are the 
Bernoulli numbers, defined as the coefficients in the Taylor expansion 



Disc 2= iLi n (z) = 27rz 



In™' 1 z 
(n-1)! ' 



(3.1) 




(3.2) 




(3.3) 



fc=0 



The functions P n (z) are known to satisfy nice functional equations, e.g., 




n > 2. 



(3.4) 
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For n = 2, eq. (3.2) reduces to the famous Bloch-Wigner dilogarithm, which satisfies in 
addition the functional equation 



P 2 (l-z) = -P 2 (z). 



(3.5) 



Just like it is possible to define single- valued analogues of the classical poly logarithms, 
there have been various attempts to define single-valued versions of multiple polyloga- 
rithms. In particular, single-valued versions have been constructed for all multiple poly- 
logarithms up to weight three [41] as well as for all harmonic polylogarithms with indices 
aj <G {0,1} [30]. The single- valued harmonic polylogarithms are characterized by symbols 
of the form (2.23), where Ui and V{ are lower weight tensors whose entries are drawn from 
the set {z, z, 1 — z, 1 — z}. Furthermore, they can always be written in the factorized form 



where the coefficients Cij are rational polynomials in multiple zeta values, and 3i and ctj 
are (possibly empty) sequences of 0's and l's. It is clear from eq. (2.23) that single- valued 
HPLs provide a natural subset of the possible functions that can appear in the e expansion 
of three-mass triangles. However, it is known that up to weight four all harmonic polylog- 
arithms with indices and 1 are expressible in terms of classical polylogarithms only [38]. 
As in this paper we are only interested in two-loop integrals for which polylogarithms of 
weight at most four can appear, we can restrict ourselves to the single-valued versions of 
the classical polylogarithms given in eq. (3.2). 

As we will see in subsequent sections, the single-valued versions of the classical poly- 
logarithms are however not sufficient to write down all three-mass triangle integrals at two 
loops. Indeed, our results show that the symbols of three mass triangles do not only have 
entries drawn from the set {z, z, 1 — z, 1 — z}, but also the entry z — z appears. We therefore 
conjecture that, to all loop orders and to all orders in e, three mass triangles are linear 
combinations of functions whose symbols 

1. have entries drawn from the set {z, z, 1 — z, 1 — z, z — z}, 

2. have a first entry which is either (zz) or (1 — z)(l — z). 

Our conjecture implies that, for each weight, we can predict a priori the space of possible 
transcendental functions that can appear in the answer. In particular, we can construct a 
basis of functions at each weight, and in practice this basis turns out to be rather small. 
The indecomposable basis elements of a given weight (i.e., basis elements which cannot be 
written as linear combinations of products of lower weights) contain a set of indecompos- 
able single- valued harmonic polylogarithms [30], augmented by some new functions whose 
symbols contain an entry equal to z — z. To our knowledge, these new functions have never 
been studied in the literature so far. In appendix A we therefore present a recursive con- 
struction for the new basis elements. As an example, the only new single-valued function 
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(3,r 3 (z),v 3 (i-z) 
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Qi(z),Qt(i-z) 


V,(z),V 4 (l- 


z),Vt{l-l/z),QZ(z) 



Table 1: Indecomposables basis elements up to weight four which can appear in the e of three 
mass triangle integrals. 



of weight three besides single-valued HPLs is 



aw = \ 
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+ - In \z\ 
4 
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1 — z 16 z 



r 2 
12 



1-z 



(3.7) 



Similar results for the new basis functions of weight four can be found in appendix A. The 
main difference between the new basis functions and the single- valued HPLs is that the new 
functions cannot be written in a factorized form (3.6), but they involve genuine multiple 
polylogarithms in (z, z). The proof that Qz{z) is indeed single- valued in the complex z plane 
follows from the construction of appendix A. We note that, as every multiple polylogarithm 
of weight at most three can be expressed through classical polylogarithms only, we could 
derive an expression for Qz{z) that does not involve any multiple polylogarithm. The 
result would be a combination of classical polylogarithms which individually have a very 
complicated branch cut structure, and the different cuts conspire such that Qz(z) is single- 
valued. We therefore prefer to present Qz(z) in the form (3.7). Furthermore note that, just 
like the singe- valued analogues of the classical polylogarithms (3.2), <2 3 (z) has a definite 
parity under complex conjugation, z <H> z. More generally, we can choose all the basis 
elements as eigenstates of the action of the Z 2 symmetry group corresponding to complex 
conjugation. The indecomposable basis elements up to weight four with given parity under 
complex conjugation are shown in tab. 1. Note that we introduce the short-hand 



V n {z) = 



2PJz) 



if n odd , 



2iP n (z) , if n even , 



(3.8) 



in order to absorb the normalization factor coming from the real and imaginary part (be- 
cause we will have to consider these functions as well in the region II, III, IV, where z and 
z are not complex conjugate to each other). 

Let us conclude this section with a discussion on how the symmetries of three-mass 
triangle integrals are implemented into the space of single- valued polylogarithms we just 



-9- 



defined. There is a natural action of the symmetric group S3 on three-mass triangles, acting 
by a permutation of the external legs. This 53 symmetry acts on the space of single- valued 
functions via 

z-tz, z-tl — z, z->l-l/z, (3.9) 

z->l/z, z ->■ 1/(1 - z) , z->z/(z-l). 

It is easy to check that the space of single- valued functions we just defined is closed under 
these transformations, because the set of irreducible polynomials {z, z, 1 — z, 1 — z, z — z} 
is invariant under the transformations (3.9) and also the first entry condition is preserved. 
As a consequence, the symmetries of a three-mass triangle graph are implemented into 
the space of functions via the functional equations arising from the transformations of the 
argument (3.9). We have worked out the relevant functional equations for all the basis 
functions up to weight four. The results are shown in appendix A. 



4. Three-mass triangles in region I 
4.1 The one- loop triangle 

In this section we compute the analytic expressions of the two-loop three-mass triangle 
master integrals in region I where X(l,u,v) < 0. The results for the other regions will be 
given in subsequent sections. 

As a warm-up, we consider the one-loop three-mass triangle and compute its Laurent 
expansion in e up to terms in the expansion of transcendental weight four. More precisely, 
we consider the integral 

TM,vl&<) = ^' J ^ k i (k+ ^ ?(k _ V2? , (4.1) 
with pf / and D = 4 — 2e. After Feynman parametrization we obtain 

Ti(pIpIpI <0 = -Or Z e)2 (-pl)' 1 ' 6 ^, 1, 1; «, v; e) , (4.2) 
where we defined 



"2, "3]u, v; e) = I I TT d Xi ^— J S 1 1 - V x t ) 
Jo Vfc=i V[Ui >) V ies J 

x (xi + x 2 + x 3 ) u ~ D (x 2 x 3 + xi x 2 u + x\ x 3 v) D l 2 ~ v , 



(4.3) 



with v = v\ + v 2 + v 3 and S any non-empty subset of {1, 2, 3} [42]. In the following we 
choose S = {1} and we obtain 



X(z/i,z/ 2 ,i/ 3 ;u,'y;e) 

(4 4) 

dx 2 dx 3 x v 2 ^ 1 x^- 1 {l + X2 + x 3 y- D (x 2 x 3 + x 2 u + x 3 v) D / 2 - / . 

'0 



I 1 

JO 
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The integral 1(1, 1, l;u,v; e) is finite as e — > and can thus be expanded in e under the 
integration sign. We perform the integration order by order in e using the method of 
ref. [43, 36, 44], i.e., we perform the Feynman parameter integrals one by one using the 
following recursive procedure: 

1. Choose a Feynman parameter Xj and write all polylogarithms in the integrand as 
multiple polylogarithms of the form G(a;xi), where a is a function of the remaining 
Feynman parameters and/or the external parameters. 

2. Use the shuffle algebra of multiple polylogarithms to replace every product of multiple 
polylogarithms in x; L by a linear combination of such functions. 

3. If all the denominators in the integral are linear, partial fraction and compute a 
primitive using the definition of multiple polylogarithms, eq. (2.10). 

4. Compute the value of the primitive at the boundaries of the integration region, and 
return to 1. 

Some comments are in order about this procedure: First, in order for the algorithm to 
converge, it is necessary to find at each step a Feynman parameter in which all the denom- 
inators are linear 2 . This condition is always satisfied in our case, as the integral (4.3) is 
two-dimensional and involves only denominators that are linear in each Feynman parame- 
ter. Second, the first step of rewriting all the polylogarithms in the integrand in the form 
G(a;Xi) involves the use functional equations among multiple polylogarithms. This task 
can easily be carried out by using the Hopf algebra structure of multiple polylogarithms, 
which allows to derive complicated functional equations among multiple polylogarithms in 
an easy way. Finally we have to address the question of how to take the limit Xi — s- oo of 
the primitive. We will illustrate this step on the example below. 

Using the algorithm we just described, we can easily obtain, at least in principle, the e 
expansion of the one-loop triangle to any order. We will illustrate this procedure in detail 
on the example of the coefficient of e° . In this case we are left with the integral 

Z(l,l > l;u,«;0)= / j— — aX2(XX3 (4.5) 

JO {1 + X 2 + X 3 )(X 2 X 3 +X2U + X 3 V) 

The integral over x 2 is trivial to perform and we get 

-7Y1 1 i n\ f°° a ln(l +x 3 ) -lnt; + ln (u + x 3 )-lnx 3 

1 1,1,1; u,v,0 = / dx 3 ■ . 4.6) 

Jo x 3 + (1 + u - v) x 3 + v 

The roots of the quadratic denominator are given by (— z, —z) defined in eq. (2.7). While in 
the present case the integration over x 3 could easily be carried out in terms of dilogarithms, 
we prefer nevertheless to compute it explicitly using the algorithm outlined above in order 
to illustrate its usage for the higher-order terms in the e expansion. Using eq. (2.12) the 
integral can be recast in the form 

rn -1 1 n ^ [°°A 9th xz)+\nu-\nv + G(-u; x 3 ) - G(0; x 3 ) 

X 1,1,1; u,v;0) = / dx 3 ■ — — . (4.7 

Jo (x 3 + z)(x 3 + z) 



2 Such integrals are called Fubini reducible in ref. [43]. 
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Note that if X(l,u,v) < 0, the integral is manifestly real and convergent for all positive 
values of u and v. Using eq. (2.10), it is easy to find a primitive of the integrand, 

p(z,z;x 3 ) = j— | - [G(-z, -u;x 3 ) - G(-z, -u;x 3 )] - lnu[G(-z;x 3 ) - G(-z;x 3 )] 
+ \nv[G(-z;x 3 )-G(-z;x 3 )] - [G(-z, -1; x 3 ) - G(-z, 0, x 3 )} 
- [G(-z,Q;x 3 ) - G(-z,-l;x 3 )}} . 

(4.8) 

Next we have to compute the limits X3 — > and x 3 — > 00 of the primitive. The first limit 
is trivial, and we obtain 

1(1, 1, 1; u, v; 0) = lim p(z, z;x 3 ). (4.9) 

In order to take the limit at infinity, we have to use the inversion relations, x 3 — > l/x 3 , for 
the multiple polylogarithms appearing inside the primitive. These identities can easily be 
derived using the Hopf algebra of multiple polylogarithms. Indeed, using the techniques 
developed in refs. [38, 39] it is easy to show that the following identity holds for x 3 > 
and z a generic complex number, 

p(z,z;x 3 ) = —^l -G (--,--; — ) +G[--, ; — ) - \nu[G(u; z) - G(u; z)] 

z — z \ z u x 3 J \ z u x 3 J 



+ G(u, 0; z) - G(u, 0; z) + In v 



f , M _I ; i - G (_i ; i 

z X 3 J \ z x 3 



-g( ) +g(--,0;—)+U 2 (z)-U 2 (z) 

V z x 3j V z x 3j 

+ G (--, -1; — ) - G ( --,0; — ) + \nu(\nz - \nz) 

\ z x 3 ) \ z x 3 ) 

— In w(ln z — In z) — — [In 2 z — In 2 z] + ln(l — z) In z — ln(l — f ) In z| . 

(4.10) 

The limit (4.9) is now trivial to take, and we obtain 

X(l,l,l;u,u;0) = — !-^< -\nu[G(u;z) - G(u; z)] + G(u, 0; z) -G(u,0;z) 



z — z 



+ Li2(z) — Li2(z) + lnn[lnz — Inz] — lnuflnz — lnz] (4.11) 
- ^ [In 2 z - In 2 z] + ln(l - z) In z - ln(l - z) In z| . 

Eq. (4.11) is the correct expression for the three-mass one-loop triangle in four dimensions. 
However, in section 3 we argued that all three-mass triangle up to two loops can be ex- 
pressed in region I through single-valued functions only. The polylogarithms appearing in 
1(1, 1, 1; n, v; 0) are manifestly odd under z z, and tab. 1 shows that there is precisely 
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one odd basis element of weight two. Indeed, if we compute the symbol of eq. (4.11) we 
get 

s[z(l,l,l;u,i,;0)] = _L-{(^)®i^|-[(l- z )(l-z)]«g,|} = J— S [2P 2 (z)] . (4.12) 

Note that we cannot add any rational multiple of Q 2 to the argument of the symbol in the 
right-hand side because of the parity of the function. We thus arrive at the conclusion that 

X(l,l,l;u,«;0) = j^V 2 (z). (4.13) 

Our strategy can immediately be extended to the higher-order terms in the Laurent 
expansion. The only technical difficulty that arises is that the x 2 integration is no longer 
trivial, but the integrand now involves powers of logarithms of the denominators. For 
example, the coefficient of e n is given by the integral 

n\ f°° dx 2 dx 3 



+ X 2 +X3)(X2X 3 +X 2 U + X 3 V) (4-14) 
x ln fc (l + x 2 + x 3 ) \n n ~ k {x2 x 3 + u x 2 + v x 3 ) . 
The logarithms can be written as multiple polylogarithms in x 2 , 



ln fc (l + x 2 +x 3 ) ln n k {x 2 x 3 + ux 2 + vx 3 ) 

= E E Q ~ ") Infe - P (! + -3) ^~ k -\v xs) 



p=0 q=Q 



(4.15) 



xG(-l-x 3 ;x 2 ) p G[-^^;x 2 

u + x 3 



The products of multiple polylogarithms can be linearized using the shuffle algebra, and we 
can thus rewrite the integrand in terms of rational functions and multiple polylogarithms 
in x 2 . If we treat x 3 as a constant, this integral can be performed in exactly the same way 
as the integral over x 3 discussed previously. 

Carrying out this procedure for the first few terms in the e expansion, we find 

T^IpIpI e)=- 2c T (-pi)" 1 "' + 2 ^ z ) 



+ e 2 



(4.16) 

(^Inulnv - ( 2 )r 2 (z) + 2 Qj(z)] + 0(e 3 ) I . 



As anticipated in section 3, the result for the one-loop three-mass triangle is expressed, at 
each order in e, as a combination of the single- valued polylogarithms shown in table 1. We 
stress that the functions appearing in eq. (4.16) form a basis of the space of functions in 
which the mass triangles naturally take values. Thus, this expression is minimal, and there 
are no further relations among these functions. 

While this is not the first time that an analytic expression for the one-loop three-mass 
triangle has been computed, we believe that our result is an improvement over existing 



- 13 - 



representations given in the literature. In refs. [28, 45] analytic expressions for the one- 
loop three-mass triangle in D = 4 dimensions in terms of dilogarithms were computed. 
However, in the computation of NNLO observables also the higher order terms in the e 
expansion were presented. In ref. [29] the one-loop three-mass triangle was expressed in 
region I, to all orders in e in terms of log-sine integrals, while in ref. [26] the same function 
was expressed up to 0(e 2 ) in region II in terms of complicated iterated integrals whose 
kernels involve inverse square roots of the Kallen function. These representations of the 
function are only valid in specific regions in (u, v) space, and the analytic continuation 
from one region to another can be very complicated. While it appears that our derivation 
makes explicit use of the fact that we work in region I, we will see in the next section that 
the analytic continuation of eq. (4.16) to all other regions is straightforward. 

Let us conclude this section by a comment about the symmetries of the one-loop 
three-mass triangle. Indeed, in section 3 we argued that the S3 symmetry permuting the 
external legs of the triangle graph are encoded into the space of functions via the functional 
equations corresponding to the transformations of the argument (3.9). It is easy to check 
using the functional equations for the basis functions, eqs. (3.5) and (3.4), as well as the 
functional equations of appendix B, that our result is invariant under z — > 1 — z and 
z — > 1/z. As these two transformations generate the whole symmetric group S3, eq. (4.16) 
is invariant under the whole action of S3, as expected. 

4.2 The two-loop master integrals 

After this warm-up, we turn in this section to the two- loop master integrals shown in fig. 1. 
The master integrals relevant for two-loop three-point off-shell functions were identified in 
ref. [26], 

d D kd D £ 1 





= e 2~f E e J 


^{plplpht) 


= e 27B£ J 






F±(pl,plpht) 





(iir D / 2 y (k + Pl ) 2 (k - P2 ) 2 1 2 {k - if ' 

d D kd D £ 1 

(iir D / 2 ) 2 (k + Pl f {k - p 2 ) 2 I 2 [{k - £) 2 } 2 ' 

d D kd D e 1 



(4.17) 



( i7r D/2)2 k 2 ( fc + pi y t 2 (fc + pi _ If (k- P3 - t) 2 ' 

d D kd D t 1 



{iir D / 2 ) 2 k 2 £ 2 (k + Pl ) 2 {l + P2 ) 2 {k - I + Pl ) 2 (l — k + p 2 ) 2 ' 
The master integral F3 was already computed in ref. [27], 

Fsiplplpl e) = 4 {-Plr l ~ 2e -^~= V, (l - £) + 0(6) . (4.18) 

Similarly, it was shown in ref. [28] that the non-planar triangle in four dimension is given 
by the square of the one-loop triangle, 



F^plplphe) = 4c 2 v pi _ A2 V 2 (zf + 0{e) . (4.19) 



( -P|)- 2 ' 2e 
(z - z) 2 
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As expected the integrals can be expressed in terms of single-valued polylogarithms. Fur- 
thermore, the symmetries of the diagrams are again entirely encoded into the functional 
equations among the basis functions. Indeed, it is easy to see from eqs. (3.4) and (3.5) 
that F4 is invariant under S3 transformations. The planar integral F3 is only invariant 
under a Z2 subgroup acting via permutation of the external legs p\ and P2, or equivalently 
z <H> 1 — z. From the inversion relation (3.4) and the fact that V±(z) is odd under z o z 
we immediately see that eq. (4.18) has the required symmetry, 

Vt (r-iy Vi (, _ _ _ P4 (_£_) _ Vi (r _ I) . (4 . 20) 

Next we turn to the computation of the two remaining master integrals F\ and F2, 
which correspond to one-loop triangles with a bubble insertion. While these integrals have 
already been considered in the literature, we believe that our results are genuinely new and 
constitute an improvement over the existing representations. Indeed, in ref. [28] analytic 
expressions for F\ and F2 were given in terms of classical polylogarithms up to weight three. 
It is however known that two-loop computations in general require polylogarithms up to 
weight four, and thus an additional order in e is required. This task was already performed 
in ref. [26] where expressions for F\ and F2 were derived in terms of complicated iterated 
integrals of weight four which are not obviously related to multiple polylogarithms. In the 
following we present analytic expressions for F\ and F2 which are given entirely in terms 
of multiple polylogarithms up to weight four. 

We start by integrating out the bubble subintegral in each of the two master integrals 
using the formula, 

J m D i 2 [-e 2 ]^ [-(k-e) 2 }^ y ' y(u 1 )t(v2)y(d-v 1 -v 2 ) 

(4.21) 

As a result, we can write each of the two master integrals as a one-loop triangle with one 
of the propagators raised to an e-dependent exponent. We obtain 

4 ^^-^h-^r-^^u^), (422) 



(2-K)e(l-25 0l «e) r(l + e) 2 T(l - e) 2 

where n £ {0,1} and l(vi, V2, ^3; u, v; e) is the Feynman-parametrized one-loop triangle 
given in eq. (4.4). As a consequence we can apply our strategy of section 4.1 to compute 
the master integrals Fi and F2. We will however deal with the two integrals separately as 
they have a different singularity structure. 

Let us start by considering the master integral F2, which corresponds to k = 1. Even 
though eq. (4.22) has an explicit pole in e, the integral 1(1 + e, 1, 1; u, v; e) is finite as e — >■ 0. 
We can thus expand under the integration sign and apply the strategy of section 4.1. We 
obtain 

1(1 + e, 1, 1; u, v; e) = f " * + V ' - e 2 (2) 1(1, 1, 1; u, v; e) + ^ P 4 ( 1 ~ ~) + 0(e 3 ) . 

(4.23) 



z \ z 
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We observe that the integral F 2 can be expressed as a linear combination of the one-loop 
triangle and the master integral F3, 



2 2 2., c r r(l + 2e), 2 . t fu~' + v~' 



F2(pi,pip 2 3 ;e) = -• 



i-Pl)- 



r(l + 2e) T(l - 2e) „ . 2 2 , , 2 , 



(4.24) 



r(i + e ) 2 r(i-e) 2 

So far we have no explanation for this decomposition of F2 into two other master integrals. 

Next we turn to the computation of the master integral F\. Unlike the previous 
case, the integral I(e, 1, 1; u, v; e) diverges as e — > 0, and so we cannot expand under 
the integration sign. Instead we perform the integral over x 3 for finite e in terms of 
hyper geometric functions, 



X(e, 1, l;u,v;e) 

,l-2e„,l-2e 



r(l-2e)r(l -e 
dx 2 < ■ 



r(2-3e) [(1 + X2 ) (v + X2) - ux 2 ] 



1-e 



?/, 



1 + X 2 



,3e-2 



(1 - 2e) (v + x 2 ) 



2 F l l,2-3e,2-2e; 



ux 2 



(1 + x 2 ) (f + x 2 ) 



(4.25) 



Note that in region I the argument of the hyper geometric function is always less than 1. 
The divergence in the integral comes entirely from the term in the first line. It is easy to 
write down a subtraction term which renders the integral finite, 



, r(l-2e)T(l-e) ,_i_ e 



and 



f 

Jo 



dx 2 Z(x 2 ; e) 



r(l-2e)r(l-e) 



(4.26) 



(4.27) 



/0 er(2-3e) 
After subtraction of the divergence, we can again expand under the integration sign and 
integrate out the Feynman parameters one-by-one. The only technical difficulty is that 
individual terms in the result have logarithmic divergences which cancel in the final answer. 
At the end of this procedure we find 



X(e, 1, l;u,v;e) 



- — < - + e(l -u-v)l(l + e,l,l;u,v) - e [2( 2 + In v In u] 

± 06 I 6 



+ e 

9 
2 



-6 V 3 (z) - 6 V 3 (1 - z) + i In u In 2 v + i In 2 u In v + 6 (3 
L \ Qt(z) + \ Qt(l -z)- ^(31nn - lnt;)^) + |(ln« - 31n V )P 3 (l - 2) ( 428 ) 



9 25 3 

- ^^M-*) 2 + — (hi 4 w + ln 4 t;) - — Inu In v (In 2 u + In 2 v) + -ln 2 uln 2 w 



- 7( 2 lnulnu - 12 Cs H uv ) + 7 U > + 0(e 4 ) . 



We see that the integral is composed of two pieces: one piece contains the master integral 
F 2 multiplied by an algebraic prefactor, while the second part is made of single-valued 
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polylogarithms that are even under z o z. This is in agreement with the result obtained 
in ref. [28] up to 0(e). Furthermore, we note that the function Q±{z) cannot be expressed 
through classical polylogarithms only, thus proving that F\ must necessarily contain mul- 
tiple polylogarithms. We refer to appendix A for a more detailed discussion of this matter. 

So far all our results are valid only in region I where A(l,u, v) < 0. In this region all 
the functions are expressed through single-valued functions in the complex variable z. In 
the next section we show how our results can be analytically continued to all the other 
regions shown in fig. 2. 

5. Analytic continuation of three-mass triangle integrals 

In the previous section we presented analytic results for all one and two-loop master inte- 
grals for three-mass triangle integrals valid in region I where A(l, u, v) < 0. In applications 
it is however important to know the expressions of the master integrals also in the other 
regions. In this section we perform the analytic continuation into the regions II, III and 
IV where X(l,u,v) > 0, and also into the boundary region V where A(l,u, v) = 0. 

5.1 Analytic continuation to the regions II, II and IV 

Let us start by by analytically continuing all the master integrals from region I where the 
Kallen function is negative to the regions II, III and IV where the Kallen function is positive 
(but non zero). The main difference with respect to region I is that z and z are real, and 
hence they are not complex conjugate to each other anymore. This implies that we will 
need a prescription of how to deal with polylogarithms that develop an imaginary part. 
Furthermore, this prescription must be such that, even though individual polylogarithms 
might develop an imaginary part, the individual basis functions V n and Q n remain real 
everywhere in the (it, v) plane. In order to derive this prescription, let us start in region I 
and then continue the function to some other region. As long as we are in region I, z and 
z are complex conjugate to each other, and without loss of generality we may assume that 

Imz > and Imz < 0. (5-1) 

If we want the analytic continuation to be smooth, we have to choose the prescription 

z — >■ z + ie and z — >■ z — ie . (5-2) 

Using this prescription we can work out all the analytic continuation formulae for the 
basis functions. The results are shown in appendix B. Note that none of the analytic 
continuation formulae shown in appendix B involves explicit factors of iir, which implies 
that all the basis functions are real everywhere in the (u,v) plane, despite the fact that 
individual polylogarithms might develop an imaginary part. This is a direct consequence 
of the single- valuedness of the basis functions in region I. 

Let us now concentrate on region II. In this region all the individual polylogarithms 
appearing inside the basis functions are real. Indeed, the basis functions have been chosen 
such that they are manifestly real if < z, z < 1. This condition is indeed satisfied, as 
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we will show in the following. It is easy to convince oneself that if one requires u, v and 
X(l,u,v) to be positive, then there are only three possible ranges for z and z, 

z,z<0, 0<z,z<l, z,z>l. (5.3) 

Since v < 1 in region II, we can see from the definition of z, eq. (2.7), that z must be 
positive in region II, and u < 1 then implies that we are in the range < z, z < 1. As the 
basis functions are manifestly real for < z, z < 1, we conclude that the expressions we 
derived for the master integrals in region I are still valid in region II. 

Let us now turn to the regions III and IV. First, we note that these regions are related 
to region I by a permutation of the external legs. In section 3 we argued that permutations 
of the external legs act on the single-valued polylogarithms via the transformations of the 
argument (3.9). For example, the mappings from region III or IV into region II are given 

by 

III < — ^11: (u,v) -> (u/v,l/v) (z,z)-*(z/(z-l),z/(z-l)), (5.4) 

IV < — >ll: (u,v) -»• (l/u,v/u) (z,z)-*(l/z,l/z). (5.5) 

The first line corresponds to the interchange of the momenta p2 and ps, while the second 
line corresponds to the interchange of the momenta p\ with p%. We can therefore always 
map the results in the regions III and IV to the region II, where we know that all the 
polylogarithms are real. This will be discussed in greater detail in the rest of this section. 
We note in passing that it is of course also possible to directly map region III to IV (and vice 
versa), and the integrals are even invariant under this transformation. This is expressed 
by interchanging u with v, or equivalently z — > 1 — z, which is just the permutation of p\ 
with p2- 

Let us now have a closer look at region III where v > 1 and z, z < 0. Indeed, if z < 0, 
we have 1 + u — v < y/X(l,u,v), and as v is positive this condition is only fulfilled for 
1 + u — v < 0. From v > u it then follows that we must be in region III. It is easy to see 
that the transformation (5.4) maps z, z < to < z/(z — l),z/(z — 1) < 1, i.e., region III 
is indeed mapped into region II. Note that the polylogarithms appearing inside the basis 
functions are also real for z, z < 0, so that the expression valid in region I is still valid in 
region III. 

Let us finally turn to region IV. As this region is related to region II via the mapping 
(z,z) — > (1 — z, 1 — z), we immediately see that in region IV we have z,z > 1. Individual 
polylogarithms appearing inside of the basis function may therefore develop imaginary 
parts, and we need to analytically continue them to region II via the transformation (5.5), 
using the prescription (5.2). The arguments of the functions V n and Q n can then be 
mapped into the unit interval by using the functional equations of appendix B. Let us 
illustrate this on the examples of the one-loop triangle T\ and the two-loop master integral 
F\. We can use the functional equations to rewrite the expressions for the integrals (4.16, 
4.28) such that z appears as 1/z in the argument of the polylogarithms. We obtain 
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1(1, 1, 1; u, v; e)| Region IV = 2- 



z — z 



V 2 (l/z) - 2e[Q 3 (l/z) + lnuV 2 (l/z) 



(5.6) 



( 2 + - In u In v ) P 2 (l /z) + 2 [- Q 4 (1/z) - 2 In n Q 3 (l /z) 
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- — ln 2 uV 2 (l/z) + -\nu\nvV 2 {l/z)\ + 0(e 6 ) } , 



X(e, 1, l;u,w;e)| Re gion iv = 1 ^ 3e + e(l - n - n)X(l + e, 1, 1; n, v)| Region iv 



+ e 



e [2( 2 + In v In n] + e 
3 



-6P 3 (l/z) -6P 3 (1 - z) + ^lnnln 2 w+ ^ In 2 nlnv + 6C3 



\ Qt(Vz) + ^ Q| (1 - z) - ^ ( 2 In 2 u + 3(3 Inn + i Inn) P 3 (l/z) 



+ -(Inn - 31nn)'P 3 (l - z) - -V 2 (l/z) 2 + — In 4 v - — lntdn 3 n - ^ Inn In 3 v 
2 2 16 8 24 

+ -ln 2 nln 2 w- C2 Inning + 3(3 Inn - 12 C3 lnv + 7C 4 l > + 0(e 4 ) . 



(5.7) 



5.2 Analytic continuation to region V 



Let us conclude this section by discussing the results for the master integral in region V 
where A(l,n, v) = 0, or equivalently z = z. While the master integrals are obviously finite 
in this region, they in general contain a prefactor ^-=, which is obviously divergent in the 
limit z — » z. This pole must therefore be spurious, and in the following we discuss how we 
can take the limit z —■ z in order to obtain the analytic expression for the master integrals 
in region V. 

As an example, consider the two-loop master integral F 3 given in eq. (4.18). In order 
to obtain the analytic expression in region V we need to take a limit of the form 

f(z,z) - f(z,z) d 



lim 

z— >z 



Z — Z 



(5.8) 



for some function f(z, z). The numerator and denominator in the left-hand-side of eq. (5.8) 
vanish simultaneously in the limit, and we can apply L'Hospital's rule to evaluate the limit 
in terms of the derivative of f(z, z). 

As an example, applying eq. (5.8) to the analytic expression for F 3 , eq. (4.18), we 
obtain, 



F 3 (z = z;e) =4 (-pi 



2\-l-2e 



6 



z(l - z) [ 
+ (z<*l-z) + 0(e) . 



Li 3 (z) - Li 2 (z) In z - \ In 2 z ln(l - z) - ^ 



2 J (5.9) 



We stress that eq. (5.9) is only valid for z < 1. The corresponding results in the other 
regions can be obtained from the inversion relations of the previous section. Furthermore, 
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note that, unlike the expression for generic (z, z), eq. (5.9) consists only of polylogarithms 
of transcendental weight three. This is an immediate consequence of the application of the 
differential operator in eq. (5.8). We checked eq. (5.9) numerically with Mathematica. 

In the same way we can obtain the expressions of all other master integrals in region V. 
In particular, all the terms that have a spurious at z = z are all proportional to either the 
one- loop triangle T\ or the two- loop ladder triangle F3 . It is therefore sufficient to compute 
the limit z — Y z of T\ in addition to the result for F3 we just obtained. The computation of 
the limit (5.8) is, however, technically more involved in this case because of the appearance 
of the new basis functions Qj,(z) and Q^(z) which contain multiple polylogarithms. Acting 
with the derivative on these multiple polylogarithms can generate new singular factors jz^, 
and we have to iterate eq. (5.8) until we reach a finite result. In the end, we arrive at 



Tl (z = z-e) = - 2c r ^LJf) (-pD" 1 "^ [ - 1 + e( Inz - 2) 

+ e 2 ( - - In 2 z + 2 In z + C2 - 4) + (z <-> 1 - z) + 0(e 3 ) 



(5.10) 



Interestingly all the polylogarithms disappear in the limit and the result is expressed 
through logarithms only. We checked the correctness of eq. (5.10) by computing the integral 
directly from the Feynman-representation (4.3) and letting u = z 2 and v = (1 — z) 2 from 
the start. We find perfect agreement with the expression (5.10) obtained by differentiating 
the result for generic (z, z). 

Having at our disposal the analytic results for T\ and F 3 in region V, we can easily 
obtain the corresponding results for all other master integrals. We start with the integral 
F4. The result is given by 



7-./- \ a 2 1 2\-2-2t( z ^ nz + (1 ~ z)ln(l - z)) 2 ( K it\ 

F 4 (z = z-e) =Ac 2 T {-p 2 ) 2 ^ z V 2(1 _ z y ^- + 0(e). (5.11) 



Next we turn to the integral F2, which is a combination of T\ and F3, 



cr £(1 + 26) 2 ( t-* + (l- z)-* 
e T(l + e) 2 { P3> V 2 



f 2 {-z = z- e) = -- ;); ; ^ (-^r ( - — '^-^ — e2 c 2 ) T 1 { Z = z-e) 

(5.12) 



r(l + 2e)T(l -2e) „ . , 9 , 

+ ' ri + «)'ri-«)' F,( " = ,; " + ° ( ^- 



Finally, the integral F\ is manifestly finite as z — > z, apart for the term proportional to F 2 . 
We can therefore simply set z = z, and express all the multiple polylogarithms in terms of 
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harmonic poly logarithms in z. We obtain 



F 1 (z = z;e) = 



+ 



z(l-z) 



2(1 - 2e) 1 - 3e 



(-pl)F 2 (z = z;e) 



4 r(l + 2e) T(l - 2e) 
2e(l-2e) T(l + e) 2 T(l - e) 2 



(-fg) 



2\-2e 



l-3e 



j-- e [C 2 + 21nzln(l-z)] 

+ e 2 [-12Li 3 (z) + 12Li 2 (z)lnz + 8 In 2 z ln(l - z) + 3 C3 



- 24Li 4 (z) + 12Li 3 (z)(31n z + ln(l - z)) - 24Li 2 (z) In 2 z 
32 

— In 3 z ln(l - z) - 12( 3 In z + In 2 z ln 2 (l - z) 



-ylnzln(l-z) + yC4 



+ (z o 1 - z) + C(e 4 ) . 



(5.13) 



As for the one-loop integral, we checked our results for F\ and F2 by directly evaluating 
the Feynman parameter integral (4.3) with u = z 2 and v = (1 — z) 2 . 



6. Conclusion 

In this paper we studied one and two-loop triangle integrals without internal masses and 
with three off-shell external legs. Based on general arguments coming from symbols and the 
Hopf algebra of multiple polylogarithms, we have shown that in region I where the Kallen 
function is negative the results can be expressed through certain single-valued multiple 
polylogarithms. Our results in this region are characterized by a strikingly simple and 
compact analytic form. In particular, it turns out that the master integral F2 can be 
written as a combination of the one-loop triangle and the two-loop master integral F3, 
at least up to the order relevant to two-loop computations. We have also studied the 
single-valued basis functions and we have shown that their functional equations almost 
trivialize the analytic continuation of our results to all other regions where the external 
masses have the same sign. We have checked that all our results satisfy the differential 
equations of ref. [26]. Furthermore, we checked that our expressions numerically 3 against 
FIESTA [46, 47] for various points in the different regions. In addition, we have compared 
our results numerically against the two-equal-mass results of ref. [50] . In all cases we found 
perfect agreement. 

While this is not the first time that these master integrals have been considered in the 
literature, we believe that our results go beyond previously known representations for these 
functions in several ways. First, while results for all master integrals expanded in e up to 
terms of transcendental weight four had already been derived in ref. [26], these expressions 
involve complicated iterated integrals which are not manifestly related to multiple poly- 
logarithms. Secondly, other representations available in the literature are either expressed 

3 All multiple polylogarithms were evaluated numerically using GiNaC [48, 49]. 
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through log-sine integrals, and thus only valid in specific kinematic regions (representa- 
tions in other regions can be obtained, but require the analytic continuation of the log-sine 
integrals) or not expanded high enough in e. All our results are expressed through multiple 
poly logarithms only and are expanded in e up to transcendental weight four. In addition, 
the functional equations among the basis functions make the analytic continuation of our 
results to other regions very simple. 

We conclude by mentioning that the approach we have used, i.e., the a priori de- 
termination the space of functions and their properties, is not restricted to single-valued 
functions and three-mass triangle integrals. Once a set of entries that can appear inside the 
symbol of a given loop integral has been determined, we can in the same way determine the 
space of functions relevant to Feynman integrals with more scales. Studies in this direction 
for two-loop box integrals with two external masses are currently under investigation. 
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A. Construction of the basis functions 

Our results for the three-mass triangles are expressed through the single-valued multiple 
polylogarithms briefly introduced in section 3. A basis for these functions is given by single- 
valued versions of HPLs, augmented by some new basis functions which involve z — z as 
an entry in its symbol. While a generic construction of single-valued HPLs is known for 
arbitrary weight [30], the new basis functions have, at least to our knowledge, never been 
studied in the literature. In this appendix we present a recursive procedure that allows to 
construct a basis for the space of single-valued functions whose symbol has entries drawn 
from the set Z = {z, z, 1 — z, 1 — z, z — z}. While we only present the construction on the 
example of this particular space of functions, we emphasize that it is easy to extend the 
construction to more general classes of functions whose symbols have prescribed entries 
and satisfy certain first entry conditions. 

Let us denote by 7~L the space spanned by multiple zeta values and polylogarithmic 
functions whose symbols have entries drawn from the set Z, and by %sv the subspace of 
H spanned by single- valued functions. As we will see below, the single- valuedness criterion 
is not entirely equivalent to the first entry condition: while single-valuedness implies the 
first entry condition, the inverse is not necessarily true. 

It is easy to check that Hsv is m fact a subalgebra of H (if we take the product of two 
functions with a prescribed set of discontinuities at most, the product will not have any 
new discontinuities). Furthermore, both % and T-Lsv are graded by the weight, i.e., they 
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can be written as direct sums of vector spaces of polylogarithms of a given weight, 

oo oo 

H = Q)H n and n S v = ^n S v,n, (A.l) 

n=0 n=0 

where we put by convention = Hsvfi = Q- We furthermore decompose each graded 
subspace into eigenspaces of complex conjugation 4 (i.e., we decompose every function into 
its real and imaginary part), 

H n = n+®n- and Hsv,n = 'Hl; Vjn ®'Hs Vjn . (A.2) 

Our goal is to find, up to a given weight (four in our case), a basis for the different graded 
subspaces. 

A first naive approach would consist in writing down the most general symbol of a 
given weight that satisfies the first entry condition and has the correct behavior under 
complex conjugation, and then to impose integrability of the symbol. Indeed, not every 
tensor corresponds to a function, but it can be shown that a tensor of the form 

^2 c l ^ii ® • • • ® u im (A. 3) 

I=(ii,...,i m ) 

is the symbol of a function if and only if the following integrability condition is fulfilled [35], 

^ ci [dkiuij A dlnu ij+1 ] ® • • • <8> ® <8> ■ ■ ■ ® u im = , (A. 4) 
i={n,...,i m ) 

where the hat indicates that the corresponding entry in the symbol is omitted. This naive 
approach however leads to two problems: 

1. For high weights, this approach is very inefficient, as it requires to impose the in- 
tegrability condition at every weight for all the components. It is more efficient to 
'recycle' information from lower weight, as dropping the last entry of an integrable 
tensor with prescribed first entry produces a tensor with the same properties, but of 
lower weight. 

2. This approach does not allow to incorporate multiple zeta values. However, as can 
be seen from the example of the single- valued HPLs of ref. [30], individual terms in 
a single-valued combination might contain non-single-valued pieces proportional to 
multiple zeta values. These terms are thus missed in this approach, and the resulting 
function will not necessarily be an element of Wsv, despite the fact that its symbol 
satisfies the first entry condition. In other words, single-valuedness implies the first 
entry condition, but the inverse is not necessarily true. A fully symbol-based approach 
cannot incorporate this delicate distinction. 

4 We could of course decompose the space of functions further into eigenspace of the S3 action (3.9). In 
the following we prefer not to do this. 
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We thus need a more refined way of constructing the basis functions. In ref. [39] it 
was shown that a way to incorporate zeta values into the symbol calculus is obtained by 
promoting the symbol map to the full coproduct of the Hopf algebra of multiple poly- 
logarithms [13]. In this setting the symbol is nothing but the maximal iteration of the 
coproduct. While T~L is obviously a Hopf algebra, its single-valued subspace T~Lsv is not, 
but we rather have 5 

Hsv ^n S v®n. (A.5) 

This follows immediately from the fact that the terms in the second factor do not need to 
satisfy the first entry condition, or more generally, from the fact that only the first factor 
in the coproduct embodies the information about the discontinuities of the functions [39]. 

Eq. (A.5) allows us to formulate a recursive procedure to construct a basis for %sv- 
The basis at weight one is obviously known, both for H\ and it single-valued subspace 
Usv,i, 

B+ = {ln\l-z\ 2 M\z\ 2 M(z-z)}, B r = {lni— Lln^}, fk , 

1-2 z (A.6) 

B + svl = {\n\l-z\ 2 M\z\ 2 }, B svl = %. 

Next, suppose that we know a basis of T-C^ Vn . We will show in the following how we can 
use this basis and recursively construct from eq. (A.5) a basis for %fy n+1 - The idea is 
simple: rather than constructing a basis for 7~i^y n , p we construct a basis of Hsv,n ® T~L\ 
instead and then 'lift' this basis to a basis of T~tgy n+1 using the integrability condition. 
Note however that in this way we can only obtain basis elements of T~i^y n+ \ which have a 
non-zero image in Hsv,n®'Hi under the coproduct. It is easy to see that those elements of 
Wsvn+i w hh vanishing image in Hsv,n ® %i are precisely the primitive elements in Hsv, 
i.e., elements £ € Hsv such that 

A(£) = l®£ + £®1. (A.7) 

We know that in the present case the primitive elements of weight n > 1 that are of interest 
are precisely the zeta values of depth one, £ n . Hence, our construction will produce a set 
of basis element to which we need to add at each weight the corresponding zeta value of 
depth one in order to get a full basis. 

Let us now assume that we are given a basis B^ Vn = {^IkkiJ 1 °f ^svn = 
dimQ'H^y n ) . Writing A n> i for the component of the coproduct that takes values in H n <g)Hi, 
and using the fact that 

K + SV,n + l ^ ^SV,n ® ^t) © ^SV,n ® «D . (A g) 

KSV,U + 1 ^ i^V,n ® WD © (KsV,n © «f ) . 

we see that B+ + U B~~ and B+~ U B~ + , with 

B^ 1,S2 = [h ® b 2 : 6i G BJV;n and fo 2 G #i 2 } , (A.9) 

5 Technically speaking, "Hsv is an "H-module. 
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form a basis for the spaces in the right-hand side of eq. (A. 8). However, not every basis 
element lies in the image of A n i. We thus need to construct the correct linear combinations. 
This is achieved by writing down the most general linear combination of vectors in each 
B'n 1,S2 and acting with A^-i^ <g) id and imposing the integrability condition (A. 4) in the 
last two factors. Formally, if we consider for example a generic linear combination £ of the 
basis elements in + U B~~ , then the requirement for £ to lie in the image of A nj i, i.e., 
for £ to be integrable, reads, 

(id®dAd)(A n _i,i ®id)(0 = 0, (A.10) 

where d denotes the usual differential on differential forms. Eq. (A. 10) gives rise to a 
linear system for the coefficients in the linear combination £, and the solution space of this 
equation is related to the image of %sy ra+1 under A n> i. The next step is simply to find 
functions whose image under A n i matches the solutions to eq. (A. 10), a step which in this 
case can easily be carried out using the methods of ref. [38]. 

The functions obtained in this way are not necessarily single- valued, despite the fact 
that their symbols satisfy the first entry condition. The reason for this lies in the fact the 
we lost track of all terms proportional to ( n when acting with A n _ lil . In the following we 
describe how we can uniquely fix this ambiguity in the present case. Indeed, as zeta values 
are real, we only have to deal with two different cases: 

1. In the odd sector, Hsvu+ii * ne terms we lost in eq. (A. 10) are necessarily of the form 

«Cnln3 + /3Cnln^4, af3e®. (A.ll) 
z 1 — z 

While these functions are not single-valued by themselves, they can appear together 
with other polylogarithms in a single-valued combination (see for example the last 
term in eq. (3.7)). It turns out that for each independent solution to eq. (A. 10) 
we can fix a and /3 in a unique way. Indeed, all the functions in T^sVn+i nave the 
property that they must vanish when z approaches the real axis, and so they must 
in particular vanish if both z and z approach the points or 1. The solutions to 
eq. (A. 10) contain terms that diverge in this limit, and a and (3 are then fixed by 
requiring the divergence to cancel. 

2. In the even sector, Hsvn+i' * ne re l evan t terms are 

aCnlnM 2 + /3Cnln|l-2| 2 + 7 C„hi(z-z), a/3 7 e<Q. (A.12) 

The first two functions are manifestly single- valued by themselves, so we do not need 
to fix the coefficients a and /3, but we simply add these functions to our basis. The 
remaining coefficient 7 is fixed in a similar way as in the odd sector, by requiring 
the function to have a smooth limit as z becomes real (note that the function is no 
longer required to vanish in this case). 

Following this procedure, we can in principle construct a basis for 7~L^gy n for every 
weight n. In practice, the combinatorics increases very quickly, which renders the procedure 



- 25 - 



quickly inefficient. The main reason for this combinatorial problem lies in the fact that 
the number of basis elements increases very quickly with n. The construction can however 
easily be carried out up to weight four, and as a result we obtain all the functions shown 
in tab. 1 (plus all possible products among these functions). The first time a single- valued 
function appears that cannot be expressed through single-valued HPLs alone is at weight 
three, where the new function Q 3 (z) appears. This function was already defined in eq. (3.7). 
At weight four, we find three new functions, two of which are simply related by replacing 
z by 1 — z. The two new independent functions read 



Qt(z) = -4 
-41n|z| 2 



: ).-.().-. \) +g(o,*,o,-,i 

z z \ z z 



z z \ z z 



G\ 0,i,0,-,l) +cfo,i,0,i,l 

z z \ z z 



+ [Li 2 (z) - Li 2 (z)] 2 + 2[Li 2 (z) + Li 2 (z)] s 



- 16[Li 4 (z) + Li 4 (z)] - 7[Li 2 (z) + Li 2 (z)] In 2 \z\ 2 + 3[Li 2 (z) - Li 2 (z)} In 
+ 3[Li 2 (z) + U 2 (z)] In |1 - z\ 2 In \z\ 2 + 6[Li 3 (l - z) + Li 3 (l - z)] In \z\ 2 



1 - z 
1 - I 



In \z\ 



1-z 



+ 18[Li 3 (z) + Li 3 (z)l In \z\ 2 - - t In 4 \z\ 2 + - t In 2 - In 2 \z\ 2 + ^ In 2 |1 - z\ 2 In 2 |z' 2 

4 4 1-z 4 



3 , 1-z 
+ - In 

2 1-z 



In II — z\ In — In \z\ 



(A.13) 



sr(z) =-g(o,o,K\,i) + cfo,o,H,ij+^ 
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Let us conclude this section by discussing whether or not the functions Q 4 (z) we just 
introduced can be expressed through classical polylogarithms alone. If we define 5 as the 
linear operator acting on tensors of rank four via 

S(a <g> b <g> c <g> d) = (a A b) A (c A d) , (A. 15) 

with aAb = a®b — b®a, then it follows from the conjecture of ref. [52] that a combination / 
of multiple polylogarithms of weight four can be expressed through classical polylogarithms 
only if and only if the symbol of / satisfies 5 [S(f)] = 0. Acting with S on Qf(z), we find 

5 [S {Qt{z))] ^0 and 5 [S {Q^(z))} = 0. (A.16) 

We thus conclude that we cannot express Q2(z) through classical polylogarithms alone. 
Q±{z) on the other hand could be expressed through classical polylogarithms only. The 
result would consist in a very complicated combination of classical polylogarithms of weight 
four, each of which has very complicated branch cuts which cancel mutually. Just like in 
the case of Q^(z), we refrain for this reason from giving the expression of Q±(z) in terms 
of classical polylogarithms. 



B. Functional equations for the basis functions 

In this section we summarize the functional equations satisfied by the single-valued func- 
tions shown in tab. 1. The functional equations we consider are those which arise through 
the action of the symmetric group S3 defined by eq. (3.9). The group S3 is generated by 
two elements, which we choose to correspond to 

z-tl/z and z -> 1 — z . (B.l) 

We therefore only show functional equations for these two transformations, and all others 
can be obtained by iterating the functional equations corresponding to transformations in 
eq. (B.l). The relevant functional equations for the single- valued analogues of the classical 
polylogarithms have already been given in eqs. (3.4) and (3.5). The functional equations 
for the new basis functions are 



Q 3 (i -z)= Qs(z) , 

03(1/*) = Q3{z)-\n\z\ 2 V 2 (z) 



(B.2) 



Q+(l/z) = Q+(z) - I In |1 - 2r| 2 In 3 |z| 2 - 12 C2 hi|l-z| 2 ln|z| 2 

-18 ln\z\ 2 P 3 (z) -20(3 In |^| 2 + | In 4 |z| 2 + 6 C2 In 2 \z\ 2 , 
QI0--Z) =QZ(z), (B.3) 
Ql(l/z) = QZ(z) - 2 In \z\ 2 Q 3 (z) + ^ In 2 \z\ 2 P 2 (z) 

+ i \n\z\ 2 In |1 - z\ 2 V 2 (z) . 
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